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ABSTRACT 

The Andromeda Galaxy Stellar Robotic Microlensing Project (The Angstrom Project) aims 
to use stellar microlensing events to trace the structure and composition of the inner regions 
of the Andromeda Galaxy (M31). We present microlensing rate and timescale predictions and 
spatial distributions for stellar and sub-stellar lens populations in combined disk and barred 
bulge models of M3 1 . We show that at least half of the stellar microlenses in and around the 
bulge are expected to have characteristic durations between 1 and 10 days, rising to as much 
as 80% for brown-dwarf dominated mass functions. These short-duration events are mostly 
missed by current microlensing surveys that are looking for Macho candidates in the M31 
dark matter halo. Our models predict that an intensive monitoring survey programme such 
as Angstrom, which will be able to detect events of durations upwards of a day, could de- 
tect around 30 events per season within ~ 5 arcminutes of the M3 1 centre, due to ordinary 
low-mass stars and remnants. This yield increases to more than 60 events for brown-dwarf 
dominated mass functions. The overall number of events and their average duration are sensi- 
tive diagnostics of the bulge mass, in particular the contribution of low-mass stars and brown 
dwarfs. The combination of an inclined disk, an offset bar-like bulge, and differences in the 
bulge and disk luminosity functions results in a four-way asymmetry in the number of events 
expected in each quadrant defined by the M31 disk axes. The asymmetry is sensitive to the 
bar prolongation, orientation and mass. 

Key words: gravitational lensing - stars: low mass, brown dwarfs - galaxies: bulges - galax- 
ies: individual: M31 - galaxies: structure 



1 INTRODUCTION 

Over the last decade the microlensing effect has been used to con- 
strain the abundance of compact dark matter in the Milky Way halo 
(e.g. lAlcock et alJl2000 A lAfonso et alJl2003al) . However, some of 
the most spectacular results have come not from the dark matter 
surveys but from stellar microlensing surveys directed towards the 
Galactic bulge, where th ousands of events involving ordinary stars 
have been detected (e.g. lAlcock et ai]|2000fJ: ludalski et alJl2003t 
lAfonso et alj2003rJ : ISumi et alj2003h* The optical depth measured 
by these surveys has shown that the Milky Way bulge is bar-like 
and oriented almost towards u s, a conclusion su pported by near- 
infrared surveys of the bulge jDwek et alJll995p . There remains 
some spread in the reported optical depth of the Milky Way bulge 
by the different survey teams, and this uncertainty partly reflects the 
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difficulty involved in separating foreground disk and bulge lens- 
ing populations, as well as the fact that the number of lines of 
sight through the bulge is severely restricted due to dust obscu- 
ration. Follow-up survey teams, which monitor events alerted in 
real-time by the main survey teams, have also obtained high signal- 
to-noise ratio lightcurves for many events, including the spectacular 
binary event EROS BLG-2000-5, which allowed an accurate mass 
determination to be made of the lens lAn et alj |2002|), as well a s 
a detailed limb-darkening profile of the source iFields et alj2003T) . 
These surve ys have also prod uced the first microlensing detection 
of a planet l Bond et al. 2004). Microlensing in our own galaxy is 
now an established tool for planetary, stellar and galactic astro- 
physics. 

With the advent of techniques like difference imag- 
ing lAlard & Luptonl Il998t> it has become possible to look 
for microlensing events in other galaxies where the sources 
are un resolved, the so-called pixe l -lensing regime. Several 
teams iPaulin-Henriksson et all l2003t ICalchi Novati et alJ l2003t 
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Figure 1. The geometry of the M31 disk and bulge. The disk of M31 is 
inclined at an angle of 77°, with the north-western side nearer to us. The 
bulge major axis is taken to be offset on the sky from the disk major axis by 
15°, as indicated by the twisting of the inner surface brightness isophotes 
IStark & Binnevll994l) . The M31 disk major and minor axes define 4 quad- 
rants, labelled Ql-4, which are referred to later in this paper. Events in- 
volving bulge lenses are expected to occur mainly in quadrants Ql and Q2 
where the disk stars lie behind the bulk of the bulge stars, whilst events in- 
volving disk lenses preferentially occur in quadrants Q3 and Q4 where they 
are in front of the bulk of the bulge stars. 



Riffeser et alJ 120031: Ide Jong et alj |2004 llMesich et all l2004 
fosh^HuT20M) "are already exploiting this and similar techniques 
to look for Macho dark matter candidates in the halo of the An- 
dromeda Galaxy (M31). In this paper we consider how an inten- 
sive stellar pixel-lensing survey of the bulge of M3 1 can be used to 
probe the stellar mass function and the inner M3 1 galactic structure. 
In some respects we have a better view of the M3 1 bulge than of our 
own. The M3 1 disk is inclined at 77° , so we have a relatively unob- 
scured view of the bulge. M31 is believed to be of earlier type than 
our own galaxy, and therefore thought to possess a more prominent 
bulge component. Intriguingly, surface brightness measurements of 
the central regions of M31 reveal a tw isting of the inner optical 
isophotes IWalterbos & Kennicuttl 19871) . consistent with a bar-like 
bulge th at is misa ligned from the disk major axis by about 15° on 
the sky fctark & Binnevll994l) . The relationship between the M31 
disk and bulge is shown schematically in FigureQ 

The Andromeda Galaxy Stellar Robotic Microlensing 
(Angstrom) Project aims to detect stellar microlensing events in 
the bulge region of M31. By obtaining upwards of three epochs per 
24-hour period the survey will be sensitive to short-duration stellar 
microlensing events due to low mass stars and brown dwarfs, which 
are missed by current M3 1 microlensing surveys looking for dark 
matter. Angstrom is using three telescopes in the Northern hemi- 
sphere at widely separated longitudes. A pilot season has already 
begun with the robotic 2m Liverpool Telescope on La Palma, the 
1.8m Doyak telescope at the Bohyunsan Observatory in Korea and 



the 2.4m Hiltner at the MDM Observatory in Arizona. This pa- 
per considers how many events a survey like Angstrom should see, 
and how their number, timescales and spatial distribution are sen- 
sitive to the low-mass end of the stellar mass functio n, and to the 
bulge mass and geometry. Previous studies (e.g. Gou ld 1 9941 iHanl 
1996) have considered the optical depth and event rate from sim- 
ple disk and spheroidal bulge populations. In this study we make 
a detailed examination of the pixel-lensing signal expected from 
a barred bulge, based on the survey sensitivity of the Angstrom 
Project. 

The format of the paper is as follows. In Section [2] we dis- 
cuss the basics of pixel-lensing theory and derive an expresion for 
the pixel-lensing rate for simple analytic distribution functions. In 
Section [3] we construct a set of five stellar microlensing models 
comprising a disk and barred bulge. The models reproduce the ob- 
served M31 surface brightness and are carefully constructed using 
consistent assumptions regarding the stellar mass function, lumi- 
nosity function and resulting mass-to-light ratio. The microlensing 
predictions from the models are presented in Section|4]and we end 
with a discussion of the results in Section|5| 



2 PIXEL-LENSING THEORY 

Pixel lensing differs from classical microlensing in that the sources 
are unresolved at baseline. Pixel-lensing events are observed as 
transient local flux enhancements against the background surface 
brightness. An event is detectable when the lens approaches suffi- 
ciently close to the line of sight to an unresolved source. When the 
lens passes within a fraction u of its Einstein radius the background 
source becomes magnified by a factor A — (u 2 + 2)/[u 2 {u 2 + 



4)] 



1/2 



which approximates to u for the high-magnification 



events that comprise the large majority of pixel-lensing lightcurves. 
A detection is possible when the excess photon count (A — l)iV* 
due to the microlensed source star exceeds the local noise, which 
for bright regions of galaxies is dominated by the galaxy count 
iVg a i, For purely Poisson noise we require 



(A — 1)N* ~ u~ 1 N, > aN^f, 



(1) 



where a = 3 is a typical theorists choice. Let us define u t to be the 
threshold impact parameter whereby events with a minimum im- 
pact parameter no < Ut satisfy equation Q and so are detectable. 
These events are detectable for a timescale 
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with 



(2) 



(3) 



the angular size of the Einstein radius, m the lens mass, d s the dis- 
tance to the source star and I the distance to the lens in units of 
d s . Here G and c denote the standard physical constants. In equa- 
tion fi is the relative proper motion of the lens across the line of 
sight and ~ #e/m is the Einstein radius crossing time. We call 
t v the visibility timescale of the pixel-lensing event. To characterise 
an event we must sample with a frequency that is much higher than 
.In this case each of our observations of the event around peak 
has a significance of at least a and the significance of the event 
detection as a whole is much greater. 

The differential rate of pixel microlensing for lenses of mass 
m against a given source star can be expressed as 
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d 5 r m = — Fd'gl 4 p 2 8Edvidpdud(j)dl, 
m 



(4) 



where F is the lens distribution function, V\ is the lens velocity 
along the line of sight, u < u t is the lens impact parameter, ip 
is the lens mass function and (f> describes the random orientation 
of the lens and source on the sky. Using equation J2j to change 
variables from p, to t v we can re-write equation J4j as 



d'Fr, 



-u t (u 2 — u 2 ) 3 ^ 2 Fd s f _E!hL \ dvidt v dudcpdl, (5) 

l \ ty J 



from which the event timescale distribution dT m /dt v can be ob- 
tained. 

For the models discussed in this paper we shall assume simple 
distribution functions comprising an analytic density function p and 
an isotropic Maxwellian velocity distribution. That is 



F 



(2n) 3 / 2 a 3 



exp 



2a 2 



(6) 



where p has a functional form yet to be specified, v t is the velocity 
across the line of sight and a is the one-dimensional velocity dis- 
persion. Equation can then be integrated to give the timescale 
distribution 



dT m 
dt v 



32 a 
m 
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K 2 -u 2 ) 1 / 2 

xI (4jS)exp(-S 2 )dldu, (7) 



where Iq(x) is the zeroth-order modified Bessel function, 7 = 



e E id s {u 2 t - ul) 1 ' 2 



/\/2at v and 5 = p sys ld a / \f2cr , with 



the 



proper motion of the line of sight. When considering a range of 
lens masses we must integrate the timescale distribution over the 
lens mass function i/>(m): 



dV f™"(dr m /dt v )4>(m)dm 



dt v 



r u </>(m) dm 



(8) 



where mi and m u are the lower and upper lens mass cut-offs, re- 
spectively. The integral of dF/dt v over t v then provides the pixel- 
lensing event rate per source, T. 

From equation Q we can see that itt depends upon both the 
source luminosity and the background surface brightness. More- 
over, dT/dt v is also sensitive to the source distance d s . Therefore 
a full evaluation of dF / dt v at a single timescale t v and along a sin- 
gle line of sight involves solving a five-dimensional integral over /, 
d s , u, m and source luminosity. An evaluation of T over the whole 
bulge region requires two further integrals over t v and galaxy sur- 
face brightness. We evaluate the integrals over I, d s , t v and surface 
brightness via direct numerical integration, whilst the integrals over 
u, m and source luminosity are performed using Monte-Carlo tech- 
niques. 

A further complication is that equation Q and subsequent for- 
mulae strictly hold only if the source star can be assumed to be 
point-like, that is to say its angular size 9, <C u8e- This approxi- 
mation may break down when bright giant stars are microlensed by 
low-mass brown dwarfs in the disk or bulge. In this finite-source 
regime there is differential magnification across the face of the 
source star and so the overall magnification is evaluated by inte- 
grating the microlensing magnification over the source. In our cal- 
culations we explicitly take account of finite source effects when 
necessary. 



3 SIMPLE MODELS FOR THE M31 INNER GALAXY 

Microlensing predictions necessarily depend upon the choice of 
model one adopts for the lens and source populations. There is a 
good deal of uncertainty in some of the key parameters, which in 
itself presents one of the prime drivers for the Angstrom Project. 
Given this uncertainty, it is sufficient for us to adopt relatively sim- 
ple models for the M3 1 stellar population that contain the key pa- 
rameters we are interested in exploring. To this end we consider 
a simple two-component model for the stellar populations in M3 1 
comprising an exponential disk and a barred bulge. 

At first sight stellar pixel lensing predictions appear to be at 
the mercy of several theoretical distributions, many of which are 
ill-constrained or unknown for M31. These include (for both the 
disk and bulge components) the density and velocity distributions, 
the stellar and remnant mass functions, and the source luminosity 
functions. However these distributions are tightly linked with one 
another and we can exploit this in order to reduce the parameter 
space. Firstly, we demand that our bulge and disk models are able to 
satisfactorily reproduce the observed M31 surface brightness pro- 
file along both the major and minor axes. Secondly, we demand 
that our assumed luminosity and mass functions are consistent with 
one another, that is to say our assumed luminosity functions fix the 
shape of the upper end of our mass functions where the bulk of the 
light is produced. Thirdly, the low mass end of our mass functions 
fix the density normalisations for the disk and bulge and thereby 
their mass-to-light ratios, since the light is fixed by the observed 
surface brightness profile. 

We proceed by investigating "heavy" and "light" disk and 
bulge models. The heavy models assume stellar mass functions that 
extend into the brown dwarf regime and that contain a significant 
density in such objects. The light models comprise just ordinary hy- 
drogen burning stars with a mass distribution similar to Milky Way 
stellar populations. All disk and bulge models also comprise a stel- 
lar remnant contribution and the disk models are assumed to com- 
prise a significant amount of gas which neither provides sources 
nor lenses for microlensing. The mass-to-light ratio implied by the 
adopted mass functions, along with the gas and remnant contri- 
butions, sets the density normalisation for each component by the 
requirement that the combined light from the disk and bulge sat- 
isfactorily reproduces the M3 1 surface brightness profile along the 
major and minor axes. 



3.1 Density and velocity distributions 

Our first component is the disk, for which we adopt a conventional 
double-exponential profile: 



Pd = Pd,o exp(-R/h) exp(—\z\/ H) 



(9) 



where pd.o is the central disk density, h = 5.8 kpc is the disk scale- 
length, H = 0.4 kpc is the scale-height, and R and z are cylindrical 
coordinates. For our heavy disk we find pa.o = 0.3 Mq pc~ 3 gives 
an i?-band mass-to-light ratio M/Lr = 3 and reproduces well the 
M31 surface brightness away from the bulge, whilst for the light 
disk model we set the normalisation a factor 3.4 times lower. 

Disk lenses and sources are assumed to orbit the disk with a 
circular speed of 235 km s _1 beyond 1.25 kpc (5 arcminutes) of 
the M31 centre. Inside 1.25 kpc the disk is taken to be a solid body 
rotator. Stars in the heavy disk are assumed to have random motions 
described by the Maxwellian velocity distribution in equation J6j, 
with a mean one-dimensional dispersion in the transverse direction 
of (t = 60 km s _1 at R = 2h, a position roughly equivalent to the 
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location of the Sun in the Milky Way; for the light disk we assume a 
dispersion that is a factor v3 lower than this. For other cylindrical 
distances R the velocity dispersion is scaled by exp(— R/2h), that 
is by the square root of the mid-plane density. 

The second component is the bulge. The twisting of the op- 
tical isophotes in the inner regions of M31 clearly indicates that 
the M31 bulge is barred and with a major axis w hich is offset by 
some 15° on the sky from the M3 1 disk major axis ( Stark & Binney 
1 1994) . as illustrated schematically in FigureQ The surface bright- 
ness profile of galactic bulges and elliptical galaxies is typi- 
cally observed to follow the de Vaucouleur exponential R 1 ^ 4 law 
( Wvse. Gil more & FranxNl997l) . However, in the inner regions a 
power-law profil e with surf ace brightness I oc R~ a is often ob- 
served I Binnev & Merrifield 1998), implying a volume density fall- 
off p oc R~( 1+a \ We shall consider two possible profiles, one an 
exponential and one a power law. Both distributions are modelled 
as bars that are offset from the M3 1 disk major axis. 

For the exponential bulge we adopt the following distribution: 

pb = pb,oexp{-[(a;b/a) 2 + {yb/qa) 2 + (z h /qa) 2 ] s }, (10) 

where a — 1 kpc is the bulge scale length, s = 0.75 is a power- 
law index, q — 0.6 defines the bar prolongation and Xb, Ub an d 
Zb are Cartesian coordinates aligned along the bar principal axes, 
with Zb normal to the disk plane. The choice of a and s provides a 
good fit to the surface brightness profile within 10 arcminutes, the 
region in which we are interested, though the fit becomes poorer 
at larger radii (see top panel of Figure|2j. For the heavy bulge we 
fix Pb,o = 12 Mq pc -3 giving M/Lr = 7 and a bulge mass of 
3 x 10 10 Mm , which is towards the upper end of favou red mass 
estimates (c.f. lKenj[l98^:IWidrow. Perrett & Suvull2003l) . For the 
light bulge the density is 3.4 times lower, which is at the lower end 
of the range of mass estimat es and is comparab le to the inferred 
mass of the Milky Way bulge iDwek et alJl995l) . 
For the power-law distribution we take 

Pb = Pb,o{l + [(x h /a) 2 + (yb/qa) 2 + (z h /qa) 2 ]}~ s/2 . (11) 

Here we adopt a = 0.75 kpc, q = 0.6 and a power-law index 
s — 3.5. The parameters s and a are chosen to allow a good match 
to the shape of the observed major- and minor-axis surface bright- 
ness profiles (bottom panel of Figure|2j.. The density normalisation 
for the power-law bulge is fixed at pb,o = 6 Mq pc~ 3 , giving a 
total mass of 1.7 x 10 10 Mq, intermediate to the heavy and light 
exponential bulge, though M/Lr = 2.1 which is the same as for 
the light exponential bulge. 

The i?-band surface brightness distribution predicted by the 
combined disk and exponential bulge model is shown in the up- 
per panel of Figure [2] whilst the lower panel shows the combined 
disk and power-law bulge profile. Also shown are i?-band surface 
brightness measurem ents along the major and mi nor axis of the 
M31 disk made by |Walterbos & Kennicutl <1987t) . The distribu- 
tions in equations J 101 and differ dramatically at large dis- 
tances from the bulge centre. The exponential profile gives a con- 
siderably poorer fit than does the power-law profile at distances 
beyond 10 arcminutes along the M31 minor axis. However, the two 
distributions are necessarily similar at smaller distances within the 
region of interest to us, with the power-law bulge having a slightly 
brighter peak. 

The bulge stars are taken to have Maxwellian random motions 
according to equation JSJ with a — 90 km s _1 for the light ex- 
ponential bulge, a factor v2 higher for the power-law bulge and 
a factor v3 higher for the heavy exponential bulge. The disper- 
sions are consistent with the range of a = 100 — 150 km s _1 
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Figure 2. Major and minor axis surface brightness profiles predicted for 
models 1—4 involving an exponential bulge (top) and for model 5 involving 
a power-law bulge (bottom). The model fits are shown by the solid lines 
whilst the crosses are i?-band data for the South- Western major axis and 
South-Eastern minor axis from Walterbos & Kennicutl 1 1987). The expo- 
nential bulge model fit is evidently poor beyond 10 arcminutes along the 
minor axis, but is a good fit within the inner bulge region of interest in this 
paper. 

observed bv lMcElro^ ll983l) . The bulge velocity dispersion could 
well differ significantly along each of the principal axes of the 
bar, so our chosen a represents the square root of the mean two- 
dimensional transverse dispersion. We al so adopt a bar pattern 
speed of 57 km s _1 kpc -1 out to 3.2 kpc IStark & Binnevlll994h . 
capping the circular speed at 180 km s _1 outside of this. 

3.2 Source luminosities 

From equation Q the intrinsic luminosity of the sources is a key 
factor in detecting pixel-lensing events. Fainter sources are more 
numerous but require higher magnification in order to be detected. 
The expected rate of pixel lensing therefore depends on the lumi- 
nosity function of the sources, but tends to be dominated by the 
brighter stars. 

The luminosity function of the M31 bulg e region has been 
explored in several studies (e.g. | Davi es 1 Frogel & Terndrud^9U 
iRich & Migtietilll995l:irablonka et alJll999UStephens et alJKOfm 
The high stellar densities cause severe blending problems, which 
hamper accurate calibrations of the bright end of the stellar lumi- 
nosity function. In so far as this can be corrected, recent studies 
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Table 1. The M31 disk-bulge models. Parameters for the density and velocity refer to the density distributions in section lTTl and the distribution 
function in equation jSJ. Parameters for the lens masses refer to the mass functions discussed in section 1331 For the disk the quoted velocity 

dispersion is normalised to its value at cylindrical distance R = 2h. Other parameters which are fixed to the same values for all models are 
discussed in sections l3~Tl and l3**3l 

Model Description Bulge parameters Disk parameters Remarks 

1 Light exponential pj,,o = 3.5 Mq pc~ 3 , <r = 90 km s^ 1 p<j q = 0-09 Mq P c_3 > <r = 35 km s _1 
bulge, light disk mi = 0.08 Mq, xi = -1.4 mi = 0.08 Mq, xi = —1.4 

a = 1 kpc, M/L R = 2.1 M/L R = 0.9 

2 Heavy exponential p^.o = 12 Mq pc -3 , a = 150 km s _1 q = 0-3 Mq pc - 3 , a = 60 km s — 1 
bulge, heavy disk m[ = 0.03 Mq, x 1 = -2.35 m\ = 0.01 Mq, x 1 = -2.35 

a = 1 kpc, M/L R = 7 M/L R = 3 

3 Light exponential pb,o = 3.5 Mq pc - 3 , a = 90 km s _1 p^ q = 0.3 Mq pc -3 , a = 60 km s — 1 
bulge, heavy disk mi = 0.08 Mq, x\ = — 1.4 mi = 0.01 Mq, x\ = -2.35 

a = 1 kpc, M/Lr = 2.1 M/L R = 3 

4 Light exponential p^.o = 4 Mq pc - 3 , a = 90 km s~ 1 p<j q = 0.3 Mq pc~ 3 , a = 60 km s — 1 Uses disk LF and 
bulge, heavy disk mi = 0.01 Mq, x\ = -2.35 mi = 0.01 Mq, x\ = -2.35 MF for both bulge 

a = 1 kpc, M /L R = 2.4 M /Lr = 3 and disk 

5 Power-law bulge, p bj0 = 6 Mq pc -3 , a = 120 km s~ x p dj0 = 0.09 Mq pc -3 , a = 35 km s _1 
light disk mi = 0.08 Mq, xi = -1.4 mi = 0.08 Mq, xi = -1.4 

a = 0.75 kpc, M/Lr = 2.1 M/Lr = 0.9 




M R —Mj m r~ m i 

Figure 3. The disk (left) and bulge (right) source star synthetic R versus 
R — I colour magnitude distributions used for our pixel-lensing calcula- 
tions. For the Angstrom Project stars above the dashed lines at Mr = 1.8 
are most likely to act as sources to pixel-lensing events as they require a 
magnification A < 30 to satisfy equation Q for a background surface 
brightness fiR = 20 mag arcsec" 2 . 

find no significant evidence for differences between the bright end 
of the st ellar lumi nosity function in the M3 1 and Milky Way bulge 
regions l Stephens et al. 2003). 

We use synthetic colour-magnitude data to model the M31 
source stars. These have been generated to match the chemical 
composition and the mass and luminosity functions observed for 
stars in the Milky Way disk and bulge and we assume they provide 
an adequate description of stars in the M31 disk and bulge. The 
synthetic datasets are based upon the theoretical stellar isochrones 
of the Padova group iGirardi & Salarij200ll) . The disk and bulge 
datasets each comprise around a million simulated stars with tabu- 
lated masses, ages, bolometric luminosities, effective temperatures 
and absolute magnitudes in several passbands. The R versus R — I 
colour-magnitude diagrams for a subset of the disk and bulge stars 
are shown in Figure|3| 

The main advantage in using synthetic datasets for microlens- 
ing calculations is that the radius of each star can be readily com- 



puted and so we can allow for finite source effects in a consistent 
way. For all but one of our microlensing models we shall assume 
the M31 bulge sources are described by the synthetic Milky Way 
bulge stars and the M31 disk sources by the synthetic Milky Way 
disk stars. For the other model we assume both bulge and disk 
sources are described by the synthetic Milky Way disk stars. This 
allows us to assess the impact of our choice of luminosity function 
upon the pixel lensing rate predictions. 

We allow for for a combined foreground and inter- 
nal M31 7?-band extinction of Ar = 0.5 mag when 
computing the mass-to-light ratio of the source stars. 
This is equivalent to three tim es the foreground value 
I Schlesel, Finkbeiner & Davis 1998) and is consistent with 
extinction measureme nts of Sb galaxies usi ng background galaxy 
counts iHolwerda, Gonzalez, Allen & van der Rruit 2005). 

3.3 Stellar and remnant mass functions 

For our light models we assume that the bulge and disk mass func- 
tions in M3 1 are similar to those of the Milky Way disk and bulge 
but do not extend into the brown dwarf regime. We assume a mass 
function of the form ip(m) = tj}» + ^ rem , where t/>„ and (/) rem 
are the mass functions of the stars and remnants, respectively. The 
remnant mass function is taken to be a superposition of three Dirac 
delta functions peaked at 0.6, 1.35 and 5 Mq, corresponding to 
white dwarfs, neutron stars and black holes, respectively. The ra- 
tio of their mass density contributions is set at 0.7:0.2:0.1, respec- 
tively. This prescription is a sligh t simplification of the remnant 
mass function advocated by Gould 1 2000). 

The stellar mass function ip* is modelled by a broken power 

law: 

_/ K(m/m t ) xi h<m<m t ) 
W ' \ K(m/m t ) X2 (m t < m < m u ) UZ; 

between lower and upper mass limits mi and m u . We set the tran- 
sition mass at m t = 0.5 Mq. The power law index X2 and m u are 
fixed to be consistent with our synthetic stellar population models. 
Accordingly we take m u = 1 Mq for bulge stars and 10 Mq for 
disk stars, whilst X2 = —2.35 for both populations. This is the case 
for both the heavy and light models. For the light bulge and disk 
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Figure 4. The pixel lensing spatial distribution for models 1-5 for events with visibility timescales t v = 1 - 100 days. The origin of each panel corresponds to 
the centre of M31 (shown by the cross) and the dotted lines running diagonally demarcate the major (upper left to lower right) and minor (lower left to upper 
right) axes of the M3 1 disk. The event rate contours are shown in black and are labelled in events per year per arcminute 2 . The green contour indicates a level 
of 1 event per year per arcminute 2 . 



we take mi = 0.08 Mq and xi = —1.4, which is c onsistent with 
obser vations of the Milky Way bulge mass function IZoccali et alJ 
2000). For the heavy models we instead let xi = X2 = —2.35 and 
extend the mass function down to mi = 0.03 Mq for the bulge and 
mi = 0.01 Mq for the disk. Different values of mi are needed in 
order to reproduce the correct M/Lr values for each component, 
though we could just as easily have used a single value of mi and 
different values for xi. 



The normalisation constant K is fixed by the overall stellar 
mass density, which contributes a fraction (1 — f lem — / gas ) of the 
overall density, where / rcm and / gas are the remnant and gas den- 
sity contributions, respectively. We assume / gas = in the bulge 
and 0.3 in the disk, whilst / rcm = 0.3 in the bulge and 0.15 in the 
disk. The values of both / gas and f Tem directly affect the computed 
M/Lr for each population as there is no light contribution from the 
remnants or the gas. They therefore affect the way in which changes 
to the stellar mass function translate into changes in the overall den- 
sity normalisation of a given component. Also, the larger the value 
of /rem then the greater the microlensing contribution from rem- 
nant lenses at the expense of stellar lenses, whereas a larger value 
of / gas reduces the overall microlensing signal since it is at the ex- 
pense of both the remnant and stellar density contributions. 



3.4 The microlensing models 

The models we use to compute the pixel-lensing rate are presented 
in Table Q Model 1 is the most conservative comprising a light 
disk and light exponential bulge. Model 2 considers a heavy disk 
and heavy exponential bulge. Model 3 examines the effect of a 
light exponential bulge combined with a heavy disk. The purpose 
of Model 4 is to allow us to explore the effect of the luminosity 
function on pixel-lensing predictions. It assumes a disk luminos- 
ity function and heavy disk mass function for both the bulge and 
disk populations, though the bulge has the same gas and remnant 
fraction as for the other bulge models. The resulting mass-to-light 
ratio is M/Lr — 2.4, slightly lower than the heavy disk value. So, 
whilst the bulge mass function in Model 4 is similar to the heavy 
bulge model, the bulge mass is closer to that of the light bulge be- 
cause of the low mass-to-light ratio of the disk-like stellar popula- 
tion. Finally, Model 5 considers a light disk and power-law bulge 
to allow a comparison between different functional forms for the 
bulge. 



4 MICROLENSING YIELDS 

We base our rate calculations on the parameters of the 
Angstrom Project, assuming a distance to M31 of 780 kpc 
(Stanek & Garnavich 1998) and a foreground extinction of Ar = 
0.17 iSchleeel. Finkbeiner & Davislfl998l) . The Angstrom Project 
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Table 2. Event yields for the Angstrom Project. (t v ) refers to the average 
visibility timescale for events in the range 1 d < t v < 100 d. N refers 
to the number of events per five-month season over an 1 1 arcminute field 
centred on the M31 bulge (as in Figure |4}- Angstrom is also using tele- 
scopes with 5 arcminute fields, for which the expected yields are about half 
that quoted below. Yields are also presented separately for the Eastern (Ql), 
Southern (Q2), Western (Q3) and Northern (Q4) quadrants defined in Fig- 
ureQ 



Model 


(M(d) 






A(l-100d) 




A r (l-10d) 






Ql 


Q2 


Q3 


Q4 


Total 




1 


16 


8 


8 


6 


5 


27 


13 


2 


7 


19 


17 


16 


12 


64 


51 


3 


14 


8 


8 


7 


5 


28 


15 


4 


7 


3 


3 


3 


3 


12 


10 


5 


13 


14 


13 


13 


10 


50 


28 



uses three telescopes, one with all arcminute field of view, the 
other two with 5 arcminute fields of view. The data is being anal - 
ysed using difference imaging techniques lAlard & Luptonlll998h . 
in which an optimal reference image is constructed from images 
taken under good seeing conditions and then subtracted from each 
target image after they have been convolved to match the seeing 
characteristics of the reference image. The resulting difference im- 
ages record only objects that have varied in their flux. For our cal- 
culations we assume observations are made in the R band and the 
reference image has a seeing of 1 arcsec. The typical .R-band ex- 
posure time is set at 600 sec and the detector has a zero point of 
R — 24.6 mag. For a Gaussian source point spread function (PSF), 
and assuming the noise within the PSF is dominated by the galactic 
background, the optimal signal-to-noise ratio is obtained by mea- 
suring the difference flux out to 71.5% of the encircled energy in the 
PSF, corresponding to a radius of 0.67 arcsec for our reference im- 
age. For guidance, the dashed lines in Figure|5|indicate the limiting 
source magnitude assuming the parameters above for the detection 
of a pixel-lensing event with A < 30 and a background surface 
brightness fiR — 20 mag arcsec -2 . 

Figure[4] shows the spatial distribution of pixel-lensing events 
with 1 d < t v < 100 d for the five models in Table Q The plots 
cover a field of view of 1 1 arcminutes centred on the bulge. The 
cross denotes the centre of M3 1 and the dotted lines show the orien- 
tation on the sky of the major and minor axes of the M31 disk. The 
contours show the expected number of events per square arcminute 
per year. For all models the event distribution is strongly dominated 
by bulge lenses, as can be seen in Figure[5] which shows the spa- 
tial distribution separately for disk and bulge lenses. The number of 
events involving disk lenses is typically around 3 events per season, 
so it is around order of magnitude lower than the bulge rate even for 
the heavy disk models. The spatial distributions in Figure |4| there- 
fore trace the bulge structure and so are noticeably rotated away 
from the disk axes for all models. Additionally, models 1-3 clearly 
show an excess of events towards the far (South-Eastern) side of the 
M3 1 disk. For these models a significant fraction of all events in- 
volve bulge lenses and disk sources, so we expect their distribution 
to be skewed towards the far disk. Models 4 and 5 are dominated 
by events involving bulge lenses and sources, so their distributions 
appear more symmetric and more closely trace the bar profile. 

Figure [6] shows the timescale distributions for models 1-5. 
The separation between the model predictions is most evident at 
timescales t v < 10 days, emphasising the need for short-timescale 
sensitivity. In general the heavier models produce a higher rate of 
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Westwards offset (arcmins) 




Westwards offset (arcmins) 

Figure 5. The pixel lensing spatial distribution in model 1 for disk (top) and 
bulge (bottom) lenses. Contours and lines are as for Figure |41 




100 



t (days) 



Figure 6. The rate-weighted pixel-lensing timescale distribution for models 
1-5 for an 1 1 arcminute field of view centred on the M31 bulge. The models 
are most distinguishable for timescales i v < 10 d. 



events and peak at shorter durations. The rate for all models tails 
off sharply for t v > 100 days. 

Table|2|shows the average duration and the number of events 
with visibility timescales between 1 and 100 d expected from a 
five-month observing season for an 1 1 arcminute field centred on 
the M31 bulge. The average duration lies between 7 and 16 days, 
with heavier bulge models giving rise to shorter events. The shorter 
durations arise because the mass function for the heavier models 



8 Kerins et al. 



extends into the low-mass brown dwarf regime, and also because 
the velocity dispersion for the heavier models is larger. 

Excluding model 4, the total number of events with durations 
between 1 and 100 d ranges from 27 up to 64, with the larger yields 
coming from the heavier, brown-dwarf rich models. The number of 
short-duration events (1 < t v < 10 d) is also shown in Table[2] 
ranging from 13 events per season for the lighter models through to 
50 events per season for the heavy exponential bulge. For the lighter 
models about half of all events have durations below 10 days, whilst 
for the heavier models they comprise up to 80% of all events. 

Model 4 assumes a disk luminosity function and heavy disk 
mass function for the bulge stars. Observations clearly exclude such 
a luminosity function for the M31 bulge, but the purpose of this 
hybrid model is simply to assess the affect of changes to the lu- 
minosity function on pixel-lensing predictions. It is evident from 
Figure|4|and Table|2|that the yield from this model is significantly 
less than for the others; it is a factor two below the conservative 
model 1. Naively this seems a very surprising result given that 
the bulge sources are brighter for model 4. Equation Q indicates 
that brighter sources can be detected at larger impact parameters 
and therefore the rate per source is larger from equation 0. How- 
ever, the low mass-to-light ratio of the disk-like stars means that 
we must adopt a brown-dwarf dominated mass function to give 
a bulge mass which is compatible with at least the lower end of 
the preferred M31 bulge mass estimates, whilst satisfying the ob- 
served M31 surface brightness. Even with a brown-dwarf domi- 
nated mass function, a heavier bulge would be too bright to be 
consistent with surface brightness observations. The combination 
of low-mass lenses and bright giant stars typically gives rise to fi- 
nite source effects which dilute the overall microlensing magnifi- 
cation. For a 0.01 Mq lens lying at a typical distance of 1 kpc in 
front of the source, a lens-source angular separation of less than 
0.04 /larcsec is required to produce a magnification A > 10. At 
the distance of M31 this corresponds to a separation of just 6 Rq, 
so finite source effects are dominant at these lens masses, resulting 
in a significant suppression of the microlensing rate for this model. 
This model is obviously rather extreme, but the basic point is that 
the boost in the microlensing rate per source due to brighter stars is 
cancelled out by the lower density normalisation required to satisfy 
surface brightness constraints. So, even if there are no constraints 
on the permitted mass of the bulge, a boost in the number of bright 
stars does not result in a gain in the overall microlensing rate. 

The number of events expected within each quadrant Ql-4 
defined in FigureQis also given in Table|2| For models 1-3 the ra- 
tio of the number of events in each quadrant is about 1.6:1.5:1.2:1 
(Q1:Q2:Q3:Q4), whilst for model 4 it is 1.2:1:1.2:1 andformodel5 
(the power-law bulge) it is 1.5:1.3:1.3:1. The number of events in 
quadrants Ql and Q2 generally outnumber those in Q3 and Q4 
because the M3 1 disk inclination skews the distribution of bulge 
lens-disk source events towards the far disk. This along with the 
offset of the barred bulge away from the disk axis introduces an 
enhanced asymmetry between events at either end of the bar in Ql 
and Q3, giving rise overall to a four- way asymmetry. The asymme- 
try is much less evident for model 4 because the source luminosity 
and lens mass functions are the same for the disk and bulge, so the 
rate of bulge lens-disk source events is similar to the rate of disk 
lens-bulge source events. For the other models, the strength of this 
four- way asymmetry depends upon the bar offset, the degree of bar 
prolongation, and the mass of the bar, so it is a good signature with 
which to probe the bar geometry. For the models 1-3, a sample of 
around 28 events would be needed in Q4 in order to detect a differ- 
ence with Ql at the 95% confidence level. Such a difference should 



therefore be detectable within three to six seasons, depending on 
the model. However, a much more efficient strategy is simply to 
undertake a likelihood analysis of the positions and timescales of 
the entire sample in order to constrain the bar parameters. 

One factor which can affect the observed event rate is extinc- 
tion. We have allowed for both foreground and internal extinction in 
calibrating our M3 1 models (see Section l3~2l . We have also allowed 
for foreground extinction in computing the M31 distance modulus, 
but w e have ignored in ternal M31 extinctio n for our rate calcula- 
tions. Holwerda, Gonzalez, Allen & van derKruitlEool have ob- 
tained internal extinction estimates for several nearby spiral galax- 
ies by studying background galaxy counts. They find that extinction 
estimates are relatively independent of the galaxy inclination angle, 
which suggests that the dust is confined to a thin sheet in the disk 
plane. In this case it is easy to quantify the affect of extinction on 
microlensing rates. The half of the source population which lies in 
front of the dust screen is unaffected by the dust. The other half is 
dimmed by a factor \Q~ 0AA R — 0.74 for an internal extinction 
of An — 0.33 mag. From equation Q we see that, for a fixed 
galaxy brightness, a 26% reduction in source flux demands a simi- 
lar reduction in the maximum impact parameter u t in order for the 
microlensed source to be detected. Since from equation {4} the mi- 
crolensing rate scales with u t then an internal M31 extinction of 
Ar = 0.33 results in an overall reduction of 13% in the observed 
rate. 



5 DISCUSSION 

The Angstrom Project is using three telescopes to conduct a high 
time resolution survey of the bulge of the Andromeda Galaxy 
(M31). The principal aim of the survey is to detect microlensing 
events with durations above one day due to low mass stars and 
brown dwarfs in the bulge of M3 1 . 

We have shown in this study that M31 presents a ripe target 
for a high time resolution microlensing survey. Using simple an- 
alytic models for the disk and bulge of M3 1 we have shown that 
a survey like Angstrom should detect around 27 events per season 
over an 1 1 arcminute field of view for a standard stellar mass func- 
tion which is truncated at the hydrogen burning limit. Half of these 
events have durations below 10 days. If the stellar mass function 
contains a substantial brown-dwarf population the expected yield 
raises to around 64 events per season, with some 80% having dura- 
tions less than 10 days. 

The underlying mass function of the lenses can be probed 
through the average duration of detected events, as well as through 
the number of events. Mass functions with dominant brown dwarf 
populations produce shorter events but also necessarily imply a 
larger mass-to-light ratio and so heavier bulge or disk mass, giv- 
ing rise to more events. 

We find that the events are dominated by bulge lenses and that 
their spatial distribution provides a good tracer of the bulge geome- 
try. In particular, if the bulge is barred and orientated away from the 
disk major axis, as suggested by surface brightness measurements, 
then this should be evident in the spatial distribution of the events. 
The importance of this is that microlensing directly traces the un- 
derlying mass distribution, rather than the light distribution, and so 
provides a completely independent probe of the underlying bar ge- 
ometry. Lastly, the combination of an offset barred bulge, a highly 
inclined disk, and different stellar luminosity functions in the disk 
and bulge is predicted to give rise to a four-way asymmetry in the 
number of events. Different numbers of events are expected in each 
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of the four quadrants defined by the major and minor axes of the 
M31 disk. 

Aside from probing the M31 bulge structure and mass func- 
tion, the high time resolution data from the Angstrom Project 
should allow for a more efficient detectio n of exotic mi crolensing 
phenomena, such as binary lensing events iBaltz & Gondolol200lh . 
than is possible with current surveys. It will also be uniquely sen- 
sitive to variable stars and transients which undergo outbursts on 
timescales of less than 10 days. 
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